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Abstract 

Following the ideas proposed in [10] and [11] on active exterior cloaking, we present 
here a systematic integral equation method to generate suitable quasistatic fields for 
cloaking, illusions and energy focusing (with given accuracy) in multiple regions of 
interests. In the quasistatic regime, the central issue is to design appropriate source 
functions for the Laplace equation so that the resulting solution will satisfy the required 
properties. We show the existence and non-uniqueness of solutions to the problem and 
study the physically relevant unique L 2 -minimal energy solution, we also provide some 
numerical evidences on the feasibility of the proposed approach. 

Key words. Field manipulation, quasistatics, layer potentials, integral equation, minimal energy 
solution, active exterior cloaking. 
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1 Introduction 



The technique of manipulating acoustic and electromagnetic fields in desired regions of space 
has been greatly advanced in the recent years, mainly due to its fascinating applications such 
as, cloaking, the creation of illusions, secret remote communication, focusing energy, and 
novel imaging techniques. The development can be roughly classified into two categories. 

The first type of techniques attempts to passively control the fields in the regions of 
interest by changing the material properties of the medium in certain surrounding regions 
while the second type of schemes focus on the active manipulation (active control) of fields 
with the help of specially designed sources. 

In [8] the authors presented the first rigorous discussion of the passive manipulation 
of fields in the context of quasistatics cloaking (see also [ 9], [30] and [31] where the in- 
variance to a change of variables is fully explained and the transformed material are fully 
described), and was later extended in [ ] to the general case of passive manipulation of 
fields in the finite frequency regime (see also the review [i] and references therein). These 
passive strategies are now known as "transformation optics". The similar strategy in con- 
text of acoustics was proposed in [6] (see also the review [3] and references therein). The 
idea behind transformation optics/acoustics is the invariance of the corresponding Dirichlet 
to Neumann-map (boundary measurements map) considered on some external boundary 
with respect to suitable change of variables which are identity on the respective boundary. 
This result implies that two different materials (the initial one and the one obtained after 
the change of variables is applied) occupying some region of space Q, will have the same 
boundary measurements maps on dQ and thus be equivalent from the point of view of an 
external observer. This leads to a long list of important applications, such as, cloaking, field 
concentrators ([36]) or field rotators, illusion optics, etc. (see [4], [3], [7], [1] and references 
therein), cloaking sensors while maintaining their sensing capability [46], [47], 

Recently, in an effort to improve accuracy and stability of the transformation optics/acoustics, 
various regularization of this scheme have been studied (see [1 1] and references therein, [26], 
1 25], [28], [34], [32], [33]). Positive results about generating broadband low-loss metamaterial 
response have been obtained in [24], [35] and a new, more stable regularization strategy was 
recently proposed in [27]. 

In a parallel direction, many researchers focused on other alternative field manipulation 
strategies. They can be grouped into two main categories, passive designs based either on 
artificial materials with extreme properties or on geometrical arguments, and active designs 
based on the active control of fields by only using antennas with no materials needed in the 
scheme, 

Among the alternative passive techniques proposed in the literature we could mention, 
plasmonic designs (see [1] and references therein), strategies based on anomalous resonance 
phenomena (see [20], [. 2], [21]), conformal mapping techniques (see [18], [17]), and compli- 
mentary media strategies (see [16]). 

Regarding the active designs for the manipulation of fields we mention that this idea 
appeared first in the context of low-frequency acoustics where various techniques for the 
active control of low-frequency sound (or active noise cancellation) were proposed in the 
literature, and we could mention here the pioneer works of Leug [44] (feedforward control of 
sound) and Olson & May [45] (feedback control of sound). For a more detailed account of 
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very interesting recent developments of the idea in the context of acoustics we mention the 
reviews [40], [42], [43], [38] [39] and the references therein. 

In the electromagnetic regime, several active designs have been recently proposed in the 
literature and we could mention the interior active cloaking strategy proposed in [19] which 
uses active boundaries and the exterior active cloaking scheme discussed in [10], [11], [12], [{ j 
(see also [48]) which uses a discrete number of active sources (antennas) to manipulate the 
fields. The active exterior strategy for 2D quasistatics cloaking was introduced in [10], were 
based on a priori information about the incoming field, with the help of one active source 
(antenna), we constructively described how one can create an almost zero field external 
region while maintaining a very small scattering effect in the far field. The proposed strategy 
did not work for objects closed to the antennas, it cloaked large objects only when they are 
far enough to the antenna (see ['•)}) and was not adaptable for three space dimensions. The 
finite frequency case was studied in the last section of [10] and in [12] (see also [9] for a recent 
review) where three active sources (antennas) were needed to create a zero field region in 
the interior of their convex hull while creating a very small scattering effect in the far field. 
The broadband character of the proposed scheme was numerically observed in [11]. We 
mention now that, from the point of view of the possible applications the constraint that 
the antennas surround the region of interests is not desirable and one would like to find a 
solution for the active manipulation of fields by using only one active source (antenna) as 
we proposed in [10]. 

In this paper, we address the problem formulated in Question 1 for the particular case of 
the quasistatic regime and a homogeneous environment). This problem is of course ill-posed 
and this explains the multitude of possible approximate solutions proposed for it. Our aim 
is to provide a unified mathematical theory for the general problem of active manipulation 
of electromagnetic or acoustic fields, which will work in a broadband regime and regardless 
of dimension, will allow for robust computational simulations and for the approximation 
of a stable optimal energy solution and will be appropriate for the more general case of 
non-homogeneous environments. In the present work we introduce the mathematical the- 
ory and analyse the problem in the quasistatic regime (modelled by the Laplace operator) 
corresponding to a homogeneous environment. 

The paper is organized as follows. In Section 2 we formulate mathematically the problem 
of generating desired field in certain regions of space using active sources. We then study 
in Section 3 the existence of solutions of the mathematical problem and Section 4 the the 
constructive approximation of a solution with minimum energy. We provide some numerical 
simulations to support our theoretical results in Section 5. Concluding and further remarks 
are offered in Section 6. Finally, for the sake of completeness, we added the proofs for two 
technical results in the Appendix. 

2 Problem formulation 

Let D s C R d (d = 2,3) be a small neighborhood of the origin 0, and D a given smooth 
domain containing Let the regions of interest, {-DfcjjtLu be N subdomains of D (i.e. 
D k CC D, 1 < k < N) that are disjoint in the sense that D k n D k > = 0,VA; ^ k! . We 
also require that Ds be disjoint with D k , Dg n D k = 0, for all k. We denote uq a smooth 
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function on M> d \D, and by a smooth function that is harmonic in a neighborhood of D^, 
i.e. Auk = in V C M. d with Dk CC V. Then the general mathematical question that we 
want to ask in quasistatic regime is: 

Question 1. Can we design an exterior active source (antenna), modeled as a continuous 
function h(x) supported on dD$, such that the harmonic field in R d \Ds generated by fa(x), 
say u, has the property that u ~ u in M\D and u ~ Uk in for all 1 < k < N , where by 
w we mean a good approximation in the uniform convergence norm? 

This question appears naturally in many applications. For instance, if the answer to 
Question 1 is positive, then one can use the active source (antenna) on D$ to generate a 
zero field in [J^ =1 and a scattering field Uq corresponding to an arbitrary object in R d \D 
to create an illusion for an external (outside of D) observer. One can also program the 
active source (antenna) to approximate N different desired fields in each of the regions D^, 
1 < k < N while creating a zero field region in M> d \D, thus sending information to regions 
of interests without being detected by an outside observer. 

We now study Question 1 in more detail. To simplify the presentation, but without loss 
of generality, we assume that all regions involved in Question 1 are balls in R d . We denote by 
B r (x) the d-dimensional open ball that centered at x G Mr with radius r > 0. Moreover, we 
first present the case where only one region of interests is involved and then, in Remark 4.1, 
show how the general result (i.e., the case of N region of interests) follows as an immediate 
consequence. Thus, let us consider Question 1 with N — 1, D$ — Bs(0), D — Br(0) and 
Di = 5 a (x ) and (for technical reasons to be discussed later) new parameters R', and a' 
such that 

a < a, R' < R, |x | > a + 5, and R' > |x | + a. (2.1) 

A schematic illustration of the problem setting and various geometrical parameters are show 
in Fig. 1. Then, in the case when uq denotes a homogeneous quasistatic potential Question 1 
can be formulated mathematically as follows. 



Formulation A. Let < e 1 be fixed. Find a function h G C(dBs(0)) such that there 
exists v G C 2 (R d \ B s (0)) n C\R d \ B 5 (0)) solution of, 



Av = in R d \ B s (0) 
v = hon dB s (0) 

\\ v ~ u l\\c(B a (x )) < ^ 
k \\ v ~ u o\\c(R d \B B XO)) — £ 



(2.2) 



where u\ is a given function harmonic in a set containing 5 a (xo) and the norm || • \\c(x) is 
the usual uniform norm on continuous functions defined on X. 

If we subtract «o from v in Formulation A, and denote by u = v — uo, g = h — uq, we 
obtain an equivalent formulation of the original problem. 
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Figure 1: The geometrical setting of Formulation A. 



Formulation A'. Let < e <C 1 be fixed. Find a function g G C(dBs(0)) such that there 
exists u G C 2 (R d \ B s (0)) n \ B s (0)) solution of, 

Aw = in R d \ B s (0) 

u = q on dBs(O) 

y V ; (2.3) 

\\U + U -ttl||c(B„(xo)) < £ 
, lhllc(K d \B H (0)) < ^ 

Thus a solution of our problem is a function g 6 C(dB$) (resp. ft, G C(dB$) for (2.2)), such 
that there exists at least a solution for problem (2.3) (resp. (2.2)). Such a solution will 
describe the required potential to be generated at the active source (antenna) so that an 
approximation of (resp. — u ) in the region £? a (x ) with e-accuracy will be possible with 
a very small perturbation of the far field (resp. very small far field). 
Let a', R', x , 5 be as before. We introduce the following space S, 

5 = L 2 (dB a ,(^)) x L 2 (dB R , m ) (2.4) 

Then H is a Hilbert space with respect to the scalar product given by 



(v,V0b= / ¥>i(y)V>i(y)<fey + / <£2(y)^2(y)ds y (2.5) 

JdB a ,(x ) JdB Rl (0) 

for all <£> = (y?i, ^2) an d ^ = ("01) ^2) in 5. The next lemma presents two technical regularity 
results which, in order to make the paper self contained, will be proved in the Appendix. 

Lemma 2.1. Let < R\ < < R2 be three constants and yo G R d an arbitrary point. Let 
f,g G C(dB Rt (y )) and definev l G C\B R ,{y ))^C l (B R ,{y Q )) andv e G C 2 (R d \B Rt (y )) n 
C 1 (R d \B Rt (y )) to be the solutions of the following interior and exterior Dirichlet problems 
respectively, 

Avi = in Br, (y ) 
Vi = f on dB Rt (y ) 
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and 



Then we have, 



[ Au e = in W l \B R Xyv) 
v e = g on dB Rt (y ) 

f 0(1) for |x| -»oo, z/d = 2 
I o(l) /or |x| — >■ 00, if d = 3 



(i) v 



zz U P 



IC(K d \B H2 (y )) 



< 



(2.7) 



where \Bi\ denotes the volume of the unit ball Bi(yo). 

The Big O and little o notations in the radiation condition guaranteeing the uniqueness 
of the solution for the exterior problem are the standard ones. 



3 Existence of solutions 



We are now ready to present the main results. Let us introduce the integral operator, 
K : L 2 (dB s (0)) ->■ E, defined as 



Ku(x,z) = (Kiu(x.),K 2 u(z)) 



(3.8) 



for any u(x) G L 2 (dBs(0)), where 



K x u{x) 
K 2 u(z) 



(9$(x, y) 

u(y) — j- 2 — ds y , for x G dB a/ (x ) 
dB s (0) av y 

u{y) ° P y) ds y , for ZG dB R ,(0) 

dB s (0) My 



(3.9) 



where v y = - — - is the normal exterior to dBs(0) and where $(x, y) represents the funda- 



mental solution of the Laplace operator, i.e., 



$(x,y) 



r 1 1 

— In- -, for d = 2 

2n |x — y| 

for d = 3 



(3.10) 



47r |x — y I 



The next result is classical but, for the sake of completeness, we included its proof in the 
Appendix. 

Lemma 3.1. The operator K defined in (3.8) is a compact linear operator from L 2 (dB$(0)) 
to E. 
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Let us introduce further the adjoint operator of K, i.e., the operator K* : H —> L 2 (dB$) 
defined through the relation, 

(Kv,u) s = (v,K*u) L 2 (dBmh Vu G E,v G L 2 (dB s (0)) (3.11) 

where (■, ■)- is the scalar product on 5 defined in (2.5) and (•, -)l 2 (9b 4 (o)) denotes the usual 
scalar product in L 2 (dBg(0)). We check, by simple change of variables and algebraic ma- 
nipulations, that the adjoint operator K* is given by, 

K*u(x)= f u 1 (y) d ^p^ds y + [ M2 ( y )^lZl dSy (3.12) 

JdB a ,{x ) VVyi JdB R ,(0) ajJ yL 

X X 

for any u — (ui, 1*2) G H and x G OBs(0), with fc» x = , — r = — . 

|x| 

From the compactness and linearity of K as given in Lemma 3.1, we conclude that the 
adjoint operator K* is compact as well. Furthermore, let us denote by Ker(K*) the kernel 
(i.e., null space) of K*. Then we have the following result. 

Proposition 3.1. Iftp — {ipi,ip2) G Ker(K*) then ip = (0,0) in 5. 

Proof. Let ip G Ker(iC*) and define 

w(x) = [ ^(y)$(x,y)^ y + / ^ 2 (y)$(x,y)rf Sy , for x6l d (3.13) 

JdB a ,( X0 ) JdB R ,(0) 

where the integrals exist as improper integrals for x G dB a i (xq) U dBn>(0). From K*ip = 
and (3.12) we have that w satisfies the Laplace equation 

Aw = 0, in B 5 (0) 

dw (3.14) 
We then conclude that 

w = constant in Bg(0) (3.15) 

We denote this constant by L, i.e., w = L in -8,5(0). Then, because by definition w is 
harmonic in Bri(0) \ B a /(pc ), from the unique continuation principle, we conclude that 

w = LmB R ,(0)\B a ,(x ) (3.16) 

The next relations for w are in fact the classical jump conditions for the single layer potentials 
with L 2 densities (see [5] and references therein). We have, 



lim / \w{x±hu x ) -w(x)| 2 ds x = (3.17) 

lim / Hx± hu x ) - w(x)| 2 rfs x = (3.18) 

ds x = (3.19) 



lim 

h ^+° JdB a ,^ ) 



lim 

h-*+0 JdB R ,(0) 



dw dw 
2 — ( x ±^ x )-2— x ±Vi X 
<w x w x 

> — ( x ±^ x )-2— x ±V 2 x 



ds x = (3.20) 



where f x = f(x) denotes the exterior normal to dB R /(0) and _B a '(x ) respectively and all 

the integral of the normal derivatives of w exists as improper integrals. From (3.16), (3.17) 
and (3.18) we obtain that 

w = L on dBjv (0) U dB a , (x ) (3.21) 

Next note that by definition w is harmonic in S a /(x ). Then, uniqueness of the interior 
Dirichlet problem for w on i? a /(x ) and (3.21) implies 

w = L in B a >(x ) (3.22) 

From (3.16), (3.22), and the two jump relations (3.19), we obtain that 

ipx = on aB a ,(x ) (3.23) 

Equation (3.23) used in the definition of w given at (3.13), implies 

w{x)= f ^ 2 (y)$(x,y)rfs y , for x6l d (3.24) 

JdB R ,(0) 

Next, relations (3.16), (3.21), and (3.22) imply that 

w = L in B R ,(0). (3.25) 
Let us now observe that Green's theorem applied to w in Br/(0) gives 

dw 

-ds* = 0. (3.26) 



dB R ,{0) 9 u x 

On the other hand, from the interior jump condition given in (3.20) together with (3.25) we 
have that 

dw I 

— = --V> 2 a.e. on dB R ,(0) (3.27) 
From (3.26) and (3.27) we deduce 

/ M*)ds x = 0. (3.28) 

JdB R ,{0) 

Observe that (3.28) guarantees the bounded behavior of w at infinity in two dimensions 
while it is well known that w will decay to zero at infinity in three dimensions. Then, the 
classical representation result for smooth functions, which are harmonic in the exterior of a 
given smooth region and bounded at infinity (see [5]), implies 

u>(x) = ^oo + / (wiy) 9 ^^ - |^(y)$(x,y)') ds y (3.29) 
JdB R ,(0) V ° u y au y J 

for all x G R d \ B R ,(0) and for some constant which depends only on the dimension. 
Using (3.21) and (3.27) in (3.29) we obtain 

f <9$(x,y) 1 f 

w(x)=w 00 + L/ — — - — ds y + - ^ 2 (y)$(x,y)rfs y 

JdB R ,(0) ° U y 1 JdB R ,(0) 

1 / \ 

= 2 Woo for xeR d \ B R ,(0) (3.30) 
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where we used (3.13) for the last integral in the first line of (3.30). Finally, (3.25) and (3.30) 
together with the pair of jump conditions given at (3.20) imply that 

^2 = a.e. on dB R ,{0). (3.31) 

The statement of the Proposition follows from (3.23) and (3.31). □ 

Before presenting the main result of this work, let us introduce the following space of 
functions 

U = K(C(dB s (0)). 
It is clear that U is a subspace of S. Moreover we have, 
Lemma 3.2. The set U C H is dense in S. 
Proof. We first observe that the subspace U C H satisfies 

U=(U ± ) ± (3.32) 

where here and further in the proof, for a given set M C H, M and M x denote its closure 
and orthogonal complement respectively in the L 2 topology generated on H by the scalar 
product defined at (2.5). Property (3.32) is classic for subspaces in a Hilbert space (see [2]). 
On the other hand we also have that 

U 1 = Ker(K*) (3.33) 

Indeed let f = (£1,62) e U ± . Then, for all up E C(dB s (0)) we have, 

= (K<p,£) s & (<p,K*£)v {aBtm = & 

& /T£ = 0<^£eKer(iT). (3.34) 

Properties (3.32) and (3.33) imply that 

U = Kex(K*) x (3.35) 

Proposition 3.1 together with (3.35) imply the density of U in 5. □ 

We are now in the position to state and prove the main result of the paper. 

Theorem 3.2. Let a,c,a',R',R be given as in (2.1). Letv = {vi,v 2 ) £ C(-B a '( x o)) x C(R d \ 
Bri(0)) be such that v\ is harmonic in _B a '(x ) and v 2 is harmonic in Mr \ B^(0). Define 
the double layer potential T> with density (p e L 2 (dB$(0)) as, 

ZV(x) = / ^(y) 9 ^ y) ds y , for x€l rf \ B s (0) 

JdB s {0) ° v y 

Then V : L 2 (dBs{0)) — > C(R d \ B a i- a+ s(0)) is a continuous operator between L 2 (dBg(0)) 
and C(IR d \ B a i- a+ $(0)) endowed with their natural topologies. Moreover, there exists a 
sequence {v n } C C(dBg(0)) such that 

Vv n — > V\ strongly in C(S a (x )), and ,T>v n — > v 2 strongly in C(R \ B R (0)) 

with respect to the uniform topology of C(B a (x )) and C(R d \ B R (0)). 
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Proof. We first observe that v G S. Then the definition of {7 and Lemma 3.2 imply that 
there exists a sequence {t> n } C (7(9^(0)) such that 

K{v n ) — > f strongly in H (3.36) 

From the definition of the H topology and (3.36) we conclude that 

\\Kiv n - ui||x 2 (as tt /(xo)) 

(3.37) 

Observe that, by definition, K±v n (resp. K 2 v n ) is the restriction to <7B a '(x ) (resp. <7Br/(0)) 
of Pw n (resp. T>v n ) where T> was defined in the statement of the Theorem. From the 
properties of V, the hypothesis on t> 1; v 2 and the regularity results of Lemma 2.1 we conclude 
that 

\\Vv n - Wl||c(B a (x )) < Cl\\ K lV n - Ui||l2(9B ,(xo)) 

(3.38) 

H^n - V 2 \\c{R d \B R (0)) < C 2 \\K 2 V n - U2 11^2(9^^,(0)) 

where we have also used the properties of a,a',R,R' and xo stated at (2.1). Finally from 
(3.37) and (3.38) we obtain the statement of the Theorem. □ 



4 The minimal energy solution 

Theorem 3.2 implies that there exist infinitely many functions h G C(dBs(0)) (resp. g G 
C(dBs(0))) as solutions to (2.2) in Formulation A ((resp. (2.3)) in Formulation A'). Indeed, 
let < e 1 and v = (vi,v 2 ) G C(5 a ,(x )) x C(R d \ B R >(0)) such that V\ is harmonic in 
S a /(xo) and v 2 is harmonic in W 1 \ Bri(0). Then, using the regularity results of Lemma 2.1 
we observe that any function h G C(dBg(0)) satisfying 

\\Kh-v\\ s <e, (4.39) 

where the || • ||a is the natural norm induced by the inner product defined in (2.5), must be 
a solution for the problem (2.2). This together with (3.38) provides a sequence of solutions 
for problem (2.2). 

Now, we will prove how, for any desired level of accuracy e, among the solutions of (4.39), 
there exists a unique solution with minimal energy norm, i.e., with minimal L 2 (dBg(0)) 
norm. We have the following result. 

Corollary 4.1. Let < e 1 and v G 2 be given. Then there exists a unique h G 
L 2 (dBs(0)) solution of the following minimization problem, 

\\ho\\L*(dB s (o)) = „ min \\h\\Li(dB s (o)) (4.40) 
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Proof. From Proposition 3.1 and classical linear operator theory we have that the linear 
bounded operator K : L 2 (dBs(0)) — >■ 5 has a dense range. This together with the classical 
theory of minimum norm solutions based on the Tikhonov regularization implies the state- 
ment of the Corollary (see [15], Theorem 16.12). In fact the classical theory implies that 
the solution h of (4.40) belongs C(dBg(0)) and is the unique solution of 

ah a + K*Kh a = K*v, with \\Kh a - v\\ s = e, (4.41) 

as the regularization strength a goes to 0. □ 

The next result is an immediate consequence of Theorem 3.2. It proves the existence of 
a class of solutions for the problem (2.3). 

Corollary 4.2. Let uq and U\ be as in (2.2) and consider v = (u\ — u ,0). Then there 
exist infinitely many functions g G C(dBs(0)) such that T>g = u with u satisfying (2.3). 
Moreover, there exists a unique function g G C(dB$(Qi)) solution of (2.3) with minimal 
L 2 (dBs(0)) norm. 

Proof. First observe that v = {u\ — u , 0) satisfies the hypothesis of Theorem 3.2 for a' > a 
satisfying (2.1) and small enough so that U\ remains harmonic on B a i(x ). Thus we have 
that there exists a sequence {g n } C C(dBs(0)) such that we have 

\\Vg n + u Q \\c{B a {* )) ^0 

(4.42) 

WDg n \\c{9. d \B R {a)) -> 

where T> is as in Theorem 3.2. Then, for < e 1 as in (2.3), we can choose N such that 
for all n > N we will have 

\\Vg n + u \\ c{ B a{xo)) < e 

(4.43) 

||'Cfi , n||c(K d \B H (0)) < e 

This implies that, there exists an index iV such that for all n > N, functions of the form T>g n 
will be solutions of the problem (2.3). Next, by using Corollary 4.1 we obtain the existence 
of a solution g G C(dBs(0)) with minimal L 2 (dBs(0)) norm. □ 

Remark 4.1. We observe that one can easily adapt the proof of Theorem 3.2 to the general 
case stated in Question 1, i.e., the case of finitely many mutually disjoint balls of inter- 
est. Thus, following the same arguments as before, one will obtain a class of solutions for 
Question 1 in this general context. Moreover, by adapting the proof of Corollary 4-1 to the 
general case of of N disjoint domains we could obtain the existence of a minimal L 2 (dBs)- 
norm solution for the problem. 

Remark 4.2. We also mention that all the results in this paper readily extend to general 
simple connected domains with C 2 boundary but for the clarity of the exposition we chose to 
present the results only in the case of spherically shaped domains. 
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Remark 4.3. It is well-known that both the interior and exterior Dirichlet problems are 
stable with respect to boundary data. This means that small perturbations on boundary data 
h produce small perturbations in the solution u in Formulation A. This further suggests that 
the inverse problem that we consider in this work, however, is unstable. To find the source 
function h that generate desired field v, we need to invert a compact integral operator. Such 
a problem is always an ill-posed problem [ , Theorem 1.17] and this is the main reason for 
the consideration of minimal energy solution. 

5 Numerical simulations 

We now present some numerical results to demonstrate the ideas that we have developed. We 
consider both two-dimensional and three-dimensional cases. To simplify the visualization, we 
only present results with regions of interests being balls, although the numerical algorithms 
we developed can deal with regions of arbitrary shapes with boundaries regular enough. The 
scattering problem (more precisely, the integral operator K) is discretized by the Nystrom 
method, following the presentation in [5] . 

In the two-dimensional case, we consider Question 1 with N — 2, 5 — 1, D\ — £?2(xi), 
T>2 = 32(^2) and D = -815(C)). The centers of D\ and D2 are xi = (0, 12) and X2 = (10,0) 
respectively. The fields are u\ = log t^t, U2 = i^p and uo = 0. The accuracy parameter 
is e = 10 _3 (||Mi|| i 2( Z ) 1 ) + ||-U2||l 2 (d 2 ) + H M o||L 2 (aD))- We show in the left plot of Fig. 2 the 
minimal energy solution of the problem with the desired fields given as above. The source 
function, supported on the unit circle, is parameterized using the azimuth angle ip G [0, 2%). 
The two middle plots of Fig. 2 show the relative differences of the field that is generated by 
the minimal energy solution and the desired field in region D\ and D2. It is clear from the 
plot that the solution strategy works almost perfectly because the mismatch between the 
desired field and the generated field is almost very small everywhere. 

In the three-dimensional case, we observe very similar results. Given an arbitrary point 
xi = (10,0,0) we considered Question 1 with N — 1, 8 = 1, Uq = 0, u\ — t^t , D± — ^(xi) 

and D = £>i 5 (0). The accuracy parameter is again e = 10~ 3 (||ui ||l 2 (Di) + ||^o||l 2 (9D))- The 
results are shown in Fig. 3. On the left plot, we show the minimal energy source h(9, ip) 
where the unit ball is parameterized using the polar angle 9 G [— 7r/2, 7r/2] and the azimuth 
angle (p G [0, 2tt). On the middle plot, we show the difference between the generated field 
and the desired filed on c?-B 2 (xi). The right plot shows the difference between the generated 
field and the desired filed on dBi^(0). Due to the limitations of visualization, we are not 
able to show the difference inside the balls which we observe to be small. 

The numerical simulations support our proof that the strategy proposed in this work on 
generating desired fields in different regions works well. More systematic numerical studies 
of the problem under various situations, together with a stability analysis will be reported 
in [7]. 
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Figure 2: Numerical results in the two-dimensional case. From top left to bottom right: the 
minimal energy source function f(ip), the relative difference between generated and desired 
fields in the neighborhoods of region Di, D2, generated field on 3D. 

6 Concluding remarks 

The idea of manipulating quasistatic fields (or in general acoustic and electromagnetic fields) 
to generate desired scattering effects has been explored extensively in the engineering com- 
munity recently due to its practical importance. In this work we present a systematic method 
to analyze mathematically and numerically the feasibility of the active field manipulation 
strategy. In the quasistatic regime, we show that one can find source functions that are able 
to generate desired quasistatic field in multiple regions of interests, to any given accuracy. 
This enable us to use the active source to create desired illusions or energy focusing without 
being detected by observations performed outside of the domain of interests. In fact, we 
show that for any given accuracy, there are infinitely many sources that can achieve the 
same effects. The source function that has the minimal energy is probably the one that 
is physically relevant and is weakly stable. Our numerical simulations confirm that the 
strategy can indeed be realized. 

The formulation that we present is independent of the spatial dimension it provides a first 
step towards the development of field manipulation techniques in more complicated settings, 
such as in low-to-medium frequency acoustic and TE or TM electromagnetic regimes even 
though the analysis in those regimes need to be done carefully due to the change of the 
integral kernels and the presence of resonances. In addition, if the problem is posed in 
a non-homogeneous medium, with known medium property, the same formulation can be 
constructed and the same type of minimal energy solution can be obtained through the 
Euler-Lagrange equation (4.41). 

Another essential discussion is about the stability of the solution. As it is well known 



13 



Figure 3: Numerical results in the three-dimensional case. Left: plot of the minimal energy 
source function f(9,ip); Middle: plot of the relative difference field on dB 2 (xi); Right: plot 
of the difference field on dBi 5 (0). 

the problem of inverting a compact operator is highly unstable and that is why we focus on 
the most physical relevant solution, namely the unique minimum energy solution. By using 
the generalized discrepancy principle it can be shown that this solution is L 2 stable with 
respect to small errors at the antennas or in the measurements of the right hand side data. 
The L 2 stability analysis together with the associated numerical discussion for the minimal 
norm solution will be presented in [37]. 

There are many potential applications of the method as we mentioned in the Introduction. 
Formulation A with u\ = corresponds to the problem of the quasistatic active exterior 
cloaking as described in [10]. It has been shown in [10, 11] (see also [12, 9] for acoustics) 
that with a few active point sources, one can generate similar effects as what we propose 
here. This is not surprising because of the non-uniqueness nature of the problem. Indeed, 
we believe that the cases in [10, 11] are special cases of the current framework, if we are 
allowed to use continuous functions to approximate the delta function model of point sources. 
Numerically, this can be done by searching for solutions with minimal L l norm instead of 
the L 2 norm (the energy norm). The numerical techniques of I 1 minimization can then be 
employed to solve the minimization problem. 
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Appendix 

Lemma 6.1. Let < R\ < < R2 be three constants and yo G M d an arbitrary point. Let 



f,g G C(&B^(y )) and definev, G C 2 ( J B«,(y ))nC 1 ( J B^(y )) and v e G C 2 (R d \B Rt (yo))H 



C 1 (M d \i?^(yo)) to be the solutions of the following interior and exterior Dirichlet problems 
respectively, 




(6.44) 
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nd 



hen we have, 



( Av e = in R d \B Rf (y ) 
v e = g on dB Rt (y ) 

f 0(1) for |x| -»oo, ifd = 2 
v e = s 

I o(l) for |x| — > 00, if d = 3 



(6.45) 



IM|c(B Hl (y )) < 



zz m f 



< 



-Bi (i?* — i?i 

R2 + -R* 



l\\g\W{dB R ,{y )) 



where \B\\ denotes the volume of the unit ball Bi(y ). 

Proof. Without loss of generality, we assume that the three balls are centered in the origin, 
i.e., yo = 0. In this condition, form the Poisson formula we have, 



and 



^.(x) 



MX] 



1 



/(y) 



II JdB Rt {0) 



g(y) 



Rl- 


|x| 2 


R* x 


-y\ d 


|x| 2 - 


-Rl 


R* x 


-y\ d 



ds v , for Ixl < R* 



d « y , 



(6.46) 



(6.47) 



where \B X \ denotes the volume of the d- dimensional unit ball. Recall that the triangle 
inequality states, 

l x — y| > ll x l - |y||> f° r an x,y e M d 



(6.48) 



From (6.46) and (6.48) we obtain 

1 



Thus 



k x < 



\B 



l/(y)l 



II JdB Rt (0) 

R* + R\ 



R* + 


X 






X 


\d-l 



\Bi\R*(R* - Ri) d 1 Jqb r ,(o) 
From (6.47) and (6.48) similarly we obtain 

1 



\g(y)\ 



II JdB Rt (0) 



R* + 


X 




R* 


X 


— R* 


d-1 



ds v , for Ixl < iL 



\f(y)\ ds yi for l x | < Rl 



ds v , for Ixl > R* 



(6.49) 



and this implies 



\vM\ < 



-R* + 


x 




\B\R* 


x 




- R* 


d-l 



! 9B Rt (0) 

By using simple algebra the last inequality becomes 



\g(y)\ds y , for |x| > i?*. 



vJx) < 



+ 



2/L 



\Bt\R. V(W -i?^- 2 (|x| -R*y-ij j 9BrA0) 
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|^(y)|rfs y , for |x| > i?*. 



and, similarly as for (6.49) this implies 

1/1 2iL 



vJx)\ < 



+ 



\B.\R, \(R 2 - R*) d - 2 (R2 - R.Y- 1 / JdBnM 

The statement of the Lemma is implied by (6.49), and (6.50). 



g(y)\ds y , for |x| > R 2 . (6.50) 



□ 



The next result is classical but for the self contained character of the paper we choose 
to include it here. 

Lemma 6.2. The operator K defined in (3.8) is a compact linear operator from L 2 (dBs(0)) 
to E. 

Proof. Let {«„} be a bounded sequence in L 2 (dBs(0)), i.e., 

\\Un\\L2(dB s (0)) < C (6.51) 

with C independent of n. Then there exists u G L 2 (dBs(0)) such that, up to a subsequence 
still indexed by n, we have 



u n — ^ u, weekly in L 2 (dB$(0)) 
Moreover, (6.51) and (6.52) imply, 

H M llL 2 (aB a (o)) < C 



(6.52) 



(6.53) 



From the continuity of the kernels of K\ and K 2 defined at (3.9) and by using (6.52) we 
obtain, 



Kiu n (x) — >■ K\u{x), for all x G dB a >(-x ) 
K 2 u n (z) -> K 2 u(z), for all z G dB R/ (0) 

Next observe that 

<9$(x,y) 1 (x-y)-y 



\Bi\ 5|x-y| 



for y G dBs(0) and for any x^y 



where again | £?i | is the volume of the (i-dimensional unit ball and we used here that 1/, 



(6.54) 
(6.55) 

y 



on dB s (0). From (6.51), (6.52), (6.53) and (6.55) we have 
\K X (u n - u)(x)| = 



<9$(x, y) 
(u n - u)(y) — ^— — rfs y 

9B 4 (0) OV y 



< \\u n — u\\Li(dB s (Q)) 



<9$(x,y) 



L 2 (dB s (0)) 



< c 



1 



|x | — a' — 5 
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d-l 



for x G <9£v(x ) 



(6.56) 



where we used that 

dist(&B a /(x ), dB s (0)) = |x | - a - 5 > 

with |xo| — a' — 5 > following from (2.4). 

Similarly as for (6.56), from (6.51), (6.52) (6.55) we obtain, 



\K 2 {u n - w)(x)| 



(« n -w)(yj— ds y 

95,5(0) c^y 



< ||« n — w||i^(aB 4 (o)) 



^(z,y) 
<9$(z,y) 



L 2 (dB s (0)) 



<C 1 for z e 0fl#(O) (6.57) 

where we used that 

dist(0B#(O), (0)) = # - 5 > 

From (6.54), (6.56) and (6.57) and the Lebesgue dominated convergence theorem, we 
obtain 

K\u n — > K\U strongly in L 2 (dB a i(x )) 

K 2 u n ^K 2 u strongly in L 2 (dB R ,(0)) (6.58) 
Then from the definition of the operator K and from (6.58) we conclude that 

Ku n — > Ku strongly in S 
and this implies the statement of the lemma. □ 
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